Adsorption of thiophene on copper surfaces: role of the dispersion force and 

electron-corehole interaction 
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We present density functional theory calculations of the geometry, adsorption energy, electronic 
density of states and bonding picture of thiophene adsorbed on Cu(lll), Cu(110) and Cu(100). 
Standard PBE functional, PBE-D and RPBE-D for including the dispersion force, and the ionic 
final-state (IFS) approximation that accounts electron-corehole interactions were employed to model 
these interfaces. The geometry was relaxed using the all four methods and the results were compared 
with the experiments. It was found that the choosing of different methods, i.e. PBE, PBE-D, 
RPBE-D, or IFS, substantially influences the geometry of thiophene/Cu(lll), but it is less sensitive 
for the other interfaces. It was interesting that RPBE-D always gives a reasonable adsorption 
energy regardless which method was adopted to relax the structure. Moreover, results also show 
that thiophene covalently bonds to all surfaces according to the analysis of electronic structures. 
All these results suggest that although theoretically predicting the interfacial structure on (111) 
surfaces is still a complicated issue, the calculations of adsorption energy should be corrected by the 
dispersion correction and standard PBE predicts correct bonding picture for these interfaces. 

PACS numbers: 68.43.Fg, 73.20.-r, 68.35.bm, 68.49.Uv 
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I. INTRODUCTION 

Molecular nanostructures on solid surfaces, especially 
on metal surfaces, have received great attention in re- 
cent years 1 - - — . Properly modeling the interfacial proper- 
ties of molecule-metal interfaces is of great importance 
to build, understand, modulate, and utilize these molec- 
ular nanostructures 1 - - — . Some building-block molecules of 
these molecular structures have small or without electri- 
cal dipolc moments, so that the dispersion force, not in- 
cluded in conventional density functional theory (DFT), 
is a major part for the van der Walls force between 
these molecules and other nanostructures, which poses 
a challenge to theoretical description of the interfacial 
properties of these molecules on solid surfaces. Profiting 
from the development of DFT, the inclusion of it was re- 
cently realized^— and applied in many molecule-metal 
interfaces 1 ^ - —. 

Despite the dispersion correction being applied to 
many interfaces, the role of it in geometry, ad- 
sorption energy, and bonding picture, however, has 
not been well recognized so far. The disper- 
sion force in Perylene-3,4,9,10-tetracarboxylic acid 
dianhydride (PTCDA)/Ag(lll), azobenzene/Ag(lll), 
thiophene/Cu(110) and pyrazine/Cu(110), was mod- 
eled using van der Walls density functional (vdW- 
DF) 14 i 15 i 19 ~ — and semi-empirical dispersion corrections 
(DFT-D)^S ~ 13 i 24 . In terms of geometry, experimental 
data of PTCDA/Ag(lll) 25 and azobenzene/Ag(lll) 26 
are available for comparison with theoretical results. It 
was found that neither vdW-DF nor DFT-D can well re- 
produce the experiment of PTCDA/Ag(lll)22, but an 
ionic final state (IFS) approximation 2 ^* 2 ^ that considers 
the electron-corehole interactions can; and only TSi 3 ., a 



variation of DFT-D, gives a reasonable Ag-N distance 
for azobenzene/Ag(lll). Although vdW-DF shows a 
good performance for the calculation of adsorption ener- 
gies, e.g. for PTCDA/Ag(lll)22.,azobenzene/Ag(lll)!£, 
and thiophene/Cu(110)^ 1 -, the performance of DFT-D, 
a much cheaper method, is still ambiguous. Further- 
more, theoretical results are even conflict to each other 
for bonding pictures. While covalent bonding picture 
was previous identified using standard DFT (PBE func- 
tional) for both pyrazine/Cu(110}22 and C 60 /Au(lll)^, 
the vdW-DF results of C 6 o/Au(lll)2£ is consistent with 
the PBE calculation, but that of pyrazine/Cu^lO)^ 2 - is, 
however, contradictive to the PBE results that the bind- 
ing of pyrazine is solely due to the dispersion force. 

Such a complex situation calls for a systematical in- 
vestigation to make an effort towards solving this is- 
sue. Since a large amount of experimental data are 
available 3 - 1 - - — , in this paper, we report density func- 
tional theory calculations of thiophene on the Cu(lll), 
Cu(110), and Cu(100) surfaces with the consideration 
of the dispersion force and electron-corehole interac- 
tions, which were, sometimes, suggested paramount in 
molecule-metal interfaces 1 ^ -28 ! 36 . The geometry of thio- 
phene adsorbed on these three Cu surfaces was relaxed 
using standard PBE, PBE-D, RPBE-D, and the IFS ap- 
proximation. Given these structures, adsorption ener- 
gies of the three interfaces were calculated using PBE, 
PBE-D, and RPBE-D, respectively, followed by a short 
summary of structural prediction for these interfaces. 
Moreover, electronic structures, including local density 
of states (LDOS), differential charge density (DCD), and 
the real space distribution of electronic states were dis- 
cussed. Finally, we conclude that, although the geometry 
issue is even more complicated, the covalent bonding is 
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primarily responsible to the molecule-metal interaction 
and the dispersion correction is helpful to correct adsorp- 
tion energy. These results provide useful guidelines for 
future applications of dispersion corrections to molecule- 
metal interfaces. 



II. COMPUTATIONAL DETAILS 

Calculations were carried out by the general gradi- 
ent approximation (GGA) for the exchange-correlation 
potential^!, the projector augmented wave metho d 38 ' 39 , 
and a plane wave basis set as implemented in the Vi- 
enna ab-initio simulation package- '^. The energy cutoff 
for plane-wave basis was set to 400 eV for all examined 
configurations. Five layers of Cu atoms, separated by a 
20-Angsrom vacuum region, were employed to model the 
Cu surfaces. Supercells p(3x3), p(2x3) and p(3x3) were 
adopted to investigate the adsorption of thiophene on 
Cu(lOO), Cu(110) and Cu(lll) surfaces, respectively. We 
anchored the S atom of thiophene on the Cu substrates 
with respect to the underlying sites, e.g., top, bridge and 
hollow sites, and rotated the molecule by the symmetry 
operations of Cu surfaces. As a result, seven, six and 
six initial configurations were considered for Cu(100), 
Cu(110) and Cu(lll), respectively. Molecules were put 
on one side of the slab with a dipole correction applied. A 
fc-mesh of 6x6x 1, verified by a 8x8 x 1 one, was adopted 
to sample all 2-D surface Brillouin zones for both geome- 
try optimization and total energy calculation. In geome- 
try optimizations, all atoms except those for the bottom 
two Cu layers were fully relaxed until the residual force 
per atom was less than 0.02 eV/A. 

There are few popular ways for including dispersion 
forces to the conventional DFT in recent years, i.e. DFT- 
D±2ri2 and vdW-DF-ii. The former is a semi-empirical 
method that the dispersion correction is parameterized 
for each element. The latter is, however, an ab-initio one, 
in which non-local correlation effects were calculated, but 
usually in a non-self-consistent manner due to huge com- 
putational demand. For the same reason, vdW-DF was 
rarely employed in structural relaxations. A large portion 
of this work is focused on structural related properties, so 
that the former method, namely DFT-D, is more suitable 
to be adopted than the latter. Structures relaxed with 
the dispersion correction, in the form of DFT-D3ii, were 
thus presented in this paper, denoted as "PBE-D". A 
different functional, namely RPBE, was also considered 
in combination with the dispersion correction, denoted 
"RPBE-D" . 

In addition, interfacial geometry of molecule-metal 
interfaces is measured by x-ray o 25 ' 26 ' 31 " — . It creates 
a core-hole when x-rays excite a core electron from 
an atom in a thiophene. Electrons nearby are prone 
to transfer to the lowest unoccupied state of the ex- 
cited molecule, screening the core-hole. In the mean- 
time, the transferred electrons start to relax to lower 
states and eventually fill the excited core level, elim- 



inating the core-hole. The electron transfer process 
was suggested much faster than the relaxation process 
of the transferred electrons in PTCDA/Ag(lll) and 
CuPcFig/Agllll) 2 ^! 2 ^. Therefore, if the excitation 
event continuously happens and the time interval of ex- 
citations is smaller than relaxation time of the trans- 
ferred electrons, some transferred electrons lack full re- 
laxation, which results in an dynamic electron accumu- 
lation around the molecule, effectively leading to the 
molecule negatively charged. The geometry of a molecule 
in such an effectively charged state, previously denoted 
IFS2£, can be very different from its neutral ground- 
state. Such a dynamic process may need time-dependent 
DFT to accurately capture, which is highly computa- 
tional demanding for these interfaces. It was how- 
ever demonstrated that the IFS approximation with a 
proper effective charge, a much cheaper method, allows 
to well reproduce the experiments of PTCDA/Ag(lll) 
and CuPcFie/Ag^ll) 2 ^^^. We thus considered the 
structural variations of thiophene adsorbed on Cu sur- 
faces with IFS effect included. Term "PBE-CH-ne" was 
used to indicate an IFS state with a total effective elec- 
tron of n transferred to the molecule, in which n can be 
a fraction. 



III. RESULTS AND DISCUSSION 

A. Structural Properties 

The results of structural relaxations were summarized 
in Tab. |U Figure [T] shows the top and side views of 
the relaxed thiophene molecules using different meth- 
ods, in which the preferred adsorption site, Cu-S bond 
length (d), and molecular tilting angle (a) were illus- 
trated. Thiophene adsorbed on each surface was sepa- 
rately discussed in this subsection, as follows. 

1. Thiophene/Cu(lll) 

Six configurations were considered in terms of adsorp- 
tion site and molecular orientation. The top site where 
the S atom sits on the top of a Cu atom, as shown in Fig. 
[lja), was found the most energetically favored among all 
the six. The Cu-S bond length (d) of the top site was cal- 
culated 3.31 A using PBE, which is 0.69 A (0.81 A) larger 
than the measured value of 2.62±0.03A^ (2.50±0.022 2 -). 
The molecular plane is roughly parallel to the surface 
with a molecular-plane-surface angle (tilting angle, a) of 
4.8°, as illustrated in Fig. [ljb). The (111) surface was 
modeled using a 3x3 supercell, equivalent to a cover- 
age of 0.11 ML. An experiment reported a tilting angle 
a = 26±5° at 0.1 ML 3 ^, which is significantly larger than 
the calculated value of 4.8°. The dispersion correction is 
able to drastically reduce the bond length to 2. 50 A and 
increase the angle to 7.0°, see Fig. [T^c). The PBE-D cal- 
culated angle is still fairly smaller than the experiment al- 
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TABLE I. Cu-S Bond lengths and molecular tilting angles of thiophene adsorbed on Cu(lll), (110), and (100) relaxed using 
PBE, PBE-D, RPBE-D, and PBE with electron corehole interaction (PBE-CH), respectively. Theoretical values of the most 
energetically favored adsorption site were reported, while two likely sites, i.e. Top and Bridge, were shown for the (100) surface. 
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FIG. 1. (Color online) (a) Adsorption site of thiophene/Cu(lll). All the four calculation methods give the same preferred 
adsorption site of thiophene/Cu(lll), in which the S is at the top of a Cu underneath and a C-S bond is oriented to the 
(110) direction of the surface, (b)-(e) Side views of the fully relaxed structures of thiophene/Cu(lll) using standard PBE (b), 
PBE-Disperion (PBE-D) (c), RPBE-D (d), and PBE with electron core-hole interaction considered (PBE-CH) (e). The values 
in side views denote the S-Cu bond length, similarly hereinafter, (f) Adsorption site of thiophene/Cu(110). (g)-(i) Side views 
of thiophene/Cu(110) relaxed using PBE (g), PBE-D (h), PBE-CH (i). (j) Two likely adsorption sites of thiophene/Cu(100). 
The configuration that the S is at the top of a Cu is slightly more stable than the other, as suggested by all the four methods. 



though the bond length is much closer. If the PBE func- 
tional was replaced by the RPBE, the result of RPBE- 
D shows a roughly flat molecule (dnpBE-D=^-00A and 
a_RPBE-n=1.4°), as shown in Fig. QJd), rather than the 
tilted molecules in those two former calculations. So far, 
none of these sets of d and a achieves a good agreement 
with the experiment. 

Three effective charges were considered in present 
work, namely 0.5e, l.Oe, and 1.5e. There is an abrupt 
change of the bond length with respect to the amount of 
effective charges. The bond length shortens to 3.15 A , 
2.59 A , and 2.42 A with effective charges of 0.5e, l.Oe, 
and 1.5e, respectively. The tilting angle a increases from 
4.8° of the standard PBE to 13.4°, 20.8°, and 20.9°, re- 
spectively (see Fig. HJe)), which appears to nearly con- 
verge at 21°. The results of an effective charge of 1.0 
e (PBE-CH-le), a very reasonable amount of charge in 
experiments, was shown in Table HI in which the bond 
length d = 2.59 A and the tilting angle a = 20.8° are 



both within the error bars of the experiment s 31 ' 33 . 



2. Thiophene/Cu(110) 

It was found that the top site is, again, the most fa- 
vorable adsorption site for thiophene/Cu(110), in which 
the Cu-S bond length is 2.42 A(PBE value) indicating an 
covalent bond. The bond length differs only O.OlA from 
the previous reported theoretical value of 2.41 A2i. DFT 
calculations revealed that thiophene has two nearly de- 
generate adsorption configurations on Cu(100), both at 
the top site, but with different tilting angles^. Similar 
to the (100) case, two nearly degenerate configurations 
were found on (110) surface. The configuration with a 
smaller tilting angle of 11° was reported in Ref. |2l|, while 
this work presented the other one with a larger angle of 
21.6°. 

According to Tab. HI either the semi-empirical disper- 
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sion correction or the corehole induced effective charge 
does not significantly change the Cu-S bond length, 
which is distinctively different from the case of thio- 
phene/Cu(lll) where the bond length appears very sen- 
sitive to the type of corrections included. The bond 
lengths revealed by PBE-D, RPBE-D, PBE-CH-le are 
2.36 A, 2.43 A, and 2.34 A, respectively. Although the 
small change of bond length, the tilting angle decreases 
from 21.6° to 11.2° and 11.8° for PBE-D and RPBE- 
D, and increases to 28.4° for PBE-CH-le, respectively. 
The increased tilting angle was ascribed to a transition 
of electronic hybridization from sp 2 to sp 3 of the carbon 
right on the top of a Cu atom, as shown in Fig. QJi), 
similar to the case of CuPcFi6/Ag(lll)21. 

There is no structural information available from ex- 
periments for this interface. The arrangement of atoms 
on Cu(llO) shares the same feature with that of Cu(100), 
i.e., the [110] directions of both surfaces are identical, 
so that the measured values of Cu(100) provide a more 
meaningful reference being compared with the calculated 
results of Cu(110) than that of Cu(lll). We thus com- 
pared our theoretical results, e.g. bond length, with the 
experimental values of Cu(100). Each theoretical value, 
in range from 2.34 A to 2.43 A, is fairly close to the ex- 
perimental value of 2.42±0.02 A^i, which implies a good 
experiment-theory agreement on the bond length. 



3. Thiophene/Cu(100) 

Earlier theoretical efforts have been made for this in- 
terface that the top site was predicted as the most favor- 
able adsorption site 4 ^. Experiments, however, suggest 
that the bridge site is more preferred. In this work, the 
top site (Fig. QTj) left) was calculated 0.1 eV more stable 
than the bridge site (Fig. [T]j) right) for all calculations 
regardless of the corrections included. The bond length 
of the top site configuration is also insensitive to differ- 
ent corrections considered, similar to the case of thio- 
phene/Cu(110), while that of the bridge site varies con- 
siderably from 2.50 A to 3.06 A. The relaxed structures 
of the top site suggest a slightly tilted molecule with the 
tilting angle from 8.7° to 10.9° except for PBE-D where 
the molecule is nearly parallel to the surface with a angle 
of 0.7°; while that of the bridge site shows carbon atoms 
lower than the sulfur giving rise to negative tilting angles 
with the only exception for PBE-CH-le. 

Our results show that theoretically predicted and ex- 
perimental measured adsorption sites are different for 
thiophene/Cu(100). According to the bond lengths 
shown in Tab. HI Cu(100) has a stronger interaction with 
thiophene than Cu(lll). It may speculate that in the 
experimental measurements 3 ^, irradiation of x-rays, to- 
gether with thermal excitation, may results in a dissocia- 
tive attachment 4 ^ 3 of thiophene on Cu(100), leaving a C4 
carbon chain adsorbed on the surface and a single S atom 
sitting at a bridge site. A dissociative attachment for a 
similar molecule comprised of phenyl- and vinyl-groups 



TABLE II. Adsorption energies of thiophene adsorbed on 
Cu(lll), (110) and (100) surfaces. Two adsorption sites were 
considered for the (100) surface. Column "Structure" indi- 
cates which functional was adopted to relax the atomic struc- 
ture and column "Energy" shows the functional employed for 
adsorption energy calculations. 

Structure Energy (111) (110) (lOO)-Top (lOO)-Bri 



PBE 


PBE 


-0.07 


-0.32 


-0.20 


-0.12 


PBE 


PBE-D 


-0.53 


-0.94 


-0.82 


-0.77 


PBE-D 


PBE-D 


-0.89 


-1.10 


-1.12 


-1.07 


PBE 


RPBE-D 


-0.42 


-0.58 


-0.54 


-0.53 


PBE-D 


RPBE-D 


-0.44 


-0.56 


-0.51 


-0.46 


RPBE-D 


RPBE-D 


-0.45 


-0.59 


-0.58 


-0.54 


Exp. 


Exp. 






-0.63 35 





and S atoms on Cu(100) was very recently observed 4 ^. It 
thus calls for further experimental and theoretical efforts 
to confirm our speculation or solve this issue. 



B. Adsorption Energy 

Adsorption energy is another important physical prop- 
erty which reflects the stability of a system. We calcu- 
lated the adsorption energies of thiophene adsorbed on 
these three surfaces using PBE, PBE-D and RPBE-D, 
and compared them with the available experimental data, 
as summarized in Tab. [TTJ 

The experimental measurement revealed that the ad- 
sorption energy for thiophene on Cu(100) is -0.63 eV 35 . 
It can be inferred from the identical [110] directions of 
the (100) and (110) surfaces and the more closely packed 
(111) surface that the adsorption energy for thiophene 
on Cu(110) should be similar to that on Cu(100), while 
that for Cu(lll) was expected a bit smaller, e.g. roughly 
-0.4—0.5 eV. Table M shows that the PBE functional, 
as expected, underestimates the adsorption energy for 
all surfaces/sites, especially for the (111) surface. Ap- 
plying the dispersion correction to the PBE functional, 
however, leads to an overestimated adsorption energy. 
It significantly improves if replacing PBE functional by 
RPBE. The RPBE-D calculated adsorption energies of 
thiophene/Cu(lll), as shown in column "(111)" , are - 
0.42 eV, -0.44 eV, and -0.45 eV for the structures relaxed 
using PBE, PBE-D and RPBE-D, respectively, consistent 
with the estimated experimental value of -0.4~-0.5 eV. 
In terms of Cu(100), the calculated adsorption energy 
using RPBE-D varies from -0.51 eV to -0.58 eV, again, 
very close to the experimental value of -0.63 eV— . It 
thus suggests that the RPBE functional combined with 
the dispersion correction is the key for properly calcu- 
lating adsorption energy, regardless which scheme, i.e. 
PBE, PBE-D, and RPBE-D, was employed to relax the 
atomistic structure. 
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C. Prediction of Interfacial Properties 

It indicates, from the three examples above, that 
the prediction of adsorption configuration and energy of 
small aromatic molecules on noble metal surfaces are very 
complicated. We intended to summarize some guidelines 
for such issue based on the experiment-theory compari- 
son of thiophene adsorbed on Cu(lll) and Cu(100) where 
experimental data are available. 

The surface electron density of Cu(lll) is rather high 
compared with other surfaces, due to the close-packed 
atomic arrangement, which reduces the corrugation of 
surface electron density. The less corrugated surface elec- 
tron density is against the formation of covalent bonds 
where a directional (corrugated) charge density is pre- 
ferred to maximize the wavefunction overlap. The ad- 
sorption energy of molecules on (111) is, therefore, a bit 
smaller than that of other surfaces. Theoretical results of 
thiophene/Cu revealed by PBE confirm this statement. 
Nevertheless, it was too weak that the thiophene-Cu(lll) 
interaction obtained by PBE (0.07 eV), ascribed to the 
lack of the non-local correlation or the dispersion correc- 
tion. The structure relaxed using PBE is also inconsis- 
tent with the experiment, roughly 0.7 A larger for the 
Cu-S distance. 

The inclusion of the dispersion correction (DFT-D) 
leads to significantly improved adsorption energies of 
thiophene on Cu surfaces, while the structural discrep- 
ancy remains unsolved. Earlier literature demonstrated 
that the electron-corehole interaction shall be consid- 
ered for theoretically reproducing the experimental struc- 
tures of two organic molecules adsorbed on Cu(lll) and 
Ag(lll) surfaces^ 2 - 7 -^^, in which conventional DFT 
faile d 23 ' 27 and vdw-DF cannot improve 2 ^. In the case 
of thiophene/Cu(lll), it was shown, again, that the con- 
sideration of electron-corehole interaction and its induced 
screening charges is crucial to compare DFT results with 
x-ray measured values. 

Electron densities of the (110) and (100) surfaces are 
relatively low and much corrugated, which gives rise to 
easier formation of typical molecule-substrate covalent 
bonds. The PBE functional suggests a pronounced co- 
valent nature of the Cu-S bonding, although the adsorp- 
tion energy was underestimated. In the presence of such 
a typical covalent bonding, the theoretical atomic struc- 
ture, e.g. the bond length, is less sensitive to the electron- 
corehole interactions, confirmed by the glycine/Cu(100) 
cas o 45 ' 46 . Therefore, the RPBE-D is sufficient to properly 
calculate the adsorption energy and might be capable to 
suitably predict the molecular configuration on these two 
surfaces. 

In a summary to this subsection, the PBE functional 
underestimates the adsorption energy for all the three 
interfaces, especially the (111) surface; while an ad- 
ditional dispersion correction to the RPBE functional 
(RPBE-D) solves this issue. The electron-corehole in- 
teraction and its induced screening charges shall be con- 
sidered in theory-experiment comparisons on the (111) 
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FIG. 2. (Color online) Local partial density of states of thio- 
phene adsorbed on Cu(lll) (a) and Cu(110) (b) calculated 
using PBE (bottom), RPBE-D (middle) and PBE-D (up). 
Variable "d" refers to the S-Cu distance, in unit of A. 



interfaces of Ag and Cu, which was repeatedly proved 
feasibl o 25 ' 27 ' 28 ' 36 . The electron-corehole interaction can 
be omitted on (110) and (100) surfaces where further 
experiments and higher level calculations shall be con- 
ducted to solve the experiment-theory discrepancy of the 
adsorption site. 



D. Electronic Structure & Bonding Picture 

1. Electronic Density of States 

It was reported that experimentally measured elec- 
tronic structures of metal-molecule interfaces are highly 
reproducible by standard DFT calculations, e.g. the 
PTCDA/Ag(lll) interface^. The dispersion correc- 
tion, as discussed in former subsections, noticeably short- 
ens the Cu-S bond length (d) and changes the molecular 
tilting angle (a) in all the four interfaces, which, most 
likely, gives rise to an appreciable variation of the elec- 
tronic structures, e.g. shifting an electronic state up- 
wards or downwards in energy. According to Tab. HI the 
dispersion correction makes thiophene 0.7 A, the largest 
change among all the interfaces, closer to the (111) sur- 
face. In addition, although PBE identifies a weaker cou- 
pling between thiophene and Cu(lll), it does indicate 
a typical covalent Cu-S bonding for thiophene/Cu(110) 
where the change of Cu-S bond length made by the dis- 
persion correction is the smallest, i.e. 0.06 A. We thus 
inspected the role of the dispersion correction to the elec- 
tronic density of states (DOS) for thiophene adsorbed on 
Cu(lll) and Cu(110), as shown in Fig. H 

Figure [2] shows that the dispersion correction does not 
significantly alter the appearance of the DOSs, especially 
for the occupied states. The 0.7 Areduction of the bond 
length lifts up those states in the energy window from - 
2.0 eV (versus the Epermi, thereinafter) to the T^Fermi by 
at most 0.2 eV, while it pushes those higher-energy states 
closer to the T^Fermi by roughly 0.6 eV. Two appreciable 
new states were found near the Epermi, from 0.2 eV to 
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FIG. 3. (Color online) Differential charge densities of PBE 
((a) and (c)) and PBE-D ((b) and (d)) calculated thio- 
phene/Cu(lll) ((a) and (b)) and thiophene/Cu(110) ((b) and 
(d)). Yellow, gray, and brown balls represent S, C, and Cu 
atoms, respectively. Slabs were cleaved along S-Cu or C-Cu 
bonds. Rainbow colormap was adopted where the red and 
blue indicate charge accumulation and reduction, respectively. 



G 

(111) (110) 

FIG. 4. (Color online) Charge density isosurfaces (PBE re- 
sults) of the hybridized LUMO-Cu states for (111) and (110) 
surfaces. Positions of C, S, and Cu atoms were marked in 
both plots. 



0.5 eV, as denoted by the two red arrows in Fig. ^a). 
In the case of thiophene/Cu(110), however, much smaller 
changes to the bond length give rise to almost unchanged 
DOSs in the energy window from -2.0 eV to 0.5 eV. The 
two pronounced peaks sitting from 1.0 to 1.5 eV split up 
more by 0.28 eV from PBE (black) to PBE-D (red), and 
a new state fades in at around 1.2 eV, as shown by the 
red arrow in Fig. [2jb). 

No appreciable charge transfer from the metal surface 
to the molecule was found for thiophene/Cu, different 
from PTCDA and other polycyclic aromatic hydrocar- 
bons, e.g. pentacene, on Cu, Ag, Au surfaces^— ' 27 i 29 ' 42 , 
in which the lowest unoccupied molecular orbital 
(LUMO) interplays with substrate states. In those sys- 
tems, a charge transfer from the substrate to the LUMO 
occurs showing a Fermi Level pinning, i.e. LUMO- 
substrate bonding state sitting just below the Fermi 
Level 6 . In terms of thiophene/Cu, the hybridized LUMO 
locates significantly higher than the Fermi Level, i.e. 
2.1 eV and 1.2 eV for the PBE calculations of thio- 
phene adsorbed on Cu(lll) and Cu(110), suggesting that 
the highest occupied molecular orbital (HOMO) may, 
most likely, play a more important role in the molecule- 
substrate electronic hybridization for thiophene on Cu 
surfaces, consistent with the electron donor nature of 
thiophene. 



2. Charge Densities 



a volume between the S and the Cu underneath (hot col- 
ors), indicating a typical Cu-S covalent bonding. A simi- 
lar Cu-C covalent bonding was also clearly illustrated for 
thiophene/Cu(110) by Fig. Efc) and (d). 

Although a dispersive bonding picture was indi- 
cated on the basis of adsorption energy by a vdW-DF 
calculation 2 !., a "covalent-like" bonding does be found in 
C6o/Au(lll) according electronic structures 2 ^. Since the 
molecule-substrate interaction of C6o/Au(lll) is believed 
weaker than that of thiophene/Cu, it strongly implies 
that a "covalent-like" bonding shall be identified by elec- 
tronic structures of the vdW-DF results of thiophene/Cu 
interfaces. In addition, according to the DCDs shown in 
Fig. EHb)-(d), the adsorption induces a dipole moment at 
the interface, which suggests an increased work-function 
of the surface upon the adsorption of thiophene. 

Furthermore, the covalent bonding picture was also 
supported by plotting the real space distribution (charge 
density), calculated using PBE, of the hybridized LUMO- 
substrate states for both the interfaces, as shown in Fig. 
21 It explicitly shows three pipelines, one for (111) and 
two for (110), of charge density connecting S or/and C 
atoms and Cu atoms underneath, which are compelling 
evidence of the S(C)-Cu covalent bonding picture for 
thiophene adsorbed on Cu surfaces. The plotted charge 
densities confirm again that the dispersion correction 
does not substantially change the covalent bonding pic- 
ture of thiophene/Cu interfaces. 



After the dispersion correction included, the electronic 
structures of these interfaces mainly keep their original 
features in terms of the appearance of DOSs. It im- 
plies that the bonding picture between thiophene and 
Cu surfaces does not be affected by the applied disper- 
sion correction, which was elucidated as follows. Figure 
[3] shows the differential charge densities (DCD), defined 

as PDCD = PThiophene/Cu ~ PThiophene ~ PCu, of PBE 

and PBE-D relaxed structures on Cu(lll) and Cu(110), 
respectively. Charge reductions (cold colors) were found 
near both Cu and S (C) atoms, especially in Fig. [SJa) 
where the reduction below the S is the most clear among 
all the four plots. These reduced charges accumulate in 



IV. CONCLUSION 

In summary, we have investigated thiophene adsorbed 
on the (111), (110), and (100) surfaces of Cu using DFT 
and post-DFT methods. It was found that standard DFT 
(PBE or RPBE) fails to predict the interfacial structure 
for thiophene/Cu(lll), underestimates the adsorption 
energy, but reveals correct bonding picture for the all in- 
terfaces. A post-DFT method, RPBE-D, solves the issue 
of underestimated adsorption energy and a calculation 
with IFS approximation perfectly predicts the interfa- 
cial structure of thiophene/Cu(lll). However, theoreti- 
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cal prediction of interfacial structures on the Cu(llO) and 
Cu(100) surface is still an issue, which needs further ex- 
periments, e.g. neutron scattering, and theory to clarify. 
In terms of electronic structures, partial DOSs of Thio- 
phene/Cu(lll) and Thiophene/Cu(110) suggest that the 
dispersion correction does not appreciably modify the oc- 
cupied states, however, slightly perturbs the unoccupied 
states, in the present cases, shifting them downwards in 
energy. Post-DFT methods, e.g. DFT-D, do not sub- 
stantially change the results of bonding picture in these 
interfaces. 

It was suggested that adopting the RPBE-D method is 
crucial for properly calculating adsorption energy and re- 
lated properties in these molecule-metal interfaces, how- 
ever, for other properties, especially electronic properties, 
standard DFT (PBE) shall give correct results. In ad- 
dition, the IFS approximation was, again, demonstrated 
feasible for predicting x-ray-measure molecular geome- 
tries on noble metal (111) surfaces. All these results shed 



considerable light on how to face with the experiment- 
theory discrepancies in molecule-metal interfaces. 
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